n=2;
X=zeros((n+1)*(2*n+1),1);
Y=zeros((n+1)*(2*n+1),1);
S=zeros((n+1)*(2*n+1),1);
C=zeros((n+1)*(2*n+1),3);
for k = -n:n
    for l = -n/2:n/2
        fprintf(1,'( %g , %g ) ',-2*l*k/n,k);
        X((l+n/2+1)+(n+1)*(k+n))=-2*l*k/n;
        Y((l+n/2+1)+(n+1)*(k+n))=k;
        S((l+n/2+1)+(n+1)*(k+n))=25+10*sign(k);
        C((l+n/2+1)+(n+1)*(k+n),:)=[0.5*abs(k)/n 0.5*abs(k)/n 0.5*abs(k)/n];
    end
    fprintf(1,'\n');
end
size(X);
size(Y);

scatter(X,Y,S,C,'filled')
axis off